BAR4R-CONF-06/028 
SLAC-PUB-12018 



A Search for B + — ► r + u Recoiling Against B — » D°£ u^X. 

The BABAR Collaboration 
February 7, 2008 



Abstract 

We present a search for the decay B + — > t + z^ t using 288 ftr 1 of data collected at the T(AS) resonance with 
the BABAR detector at the SLAC PEP-H B -Factory. A sample of events with one reconstructed semilep tonic 
B decay (B~ — > D°l~v^X) is selected, and in the recoil a search for B + — > t + v t signal is performed. The 
r is identified in the following channels: r + — > e + u e v T , r + — > ^ + f A1 z/ T , r + — > 7r + P T and r + — > 7r + 7r°^ T . 
We measure a branching fraction of ->• r+^ r ) = (0.88lo;gf (stat.) ± 0.11 (syst.)) x 10~ 4 and extract 

an upper limit on the branching fraction, at the 90% confidence level, of B{B + — > t + v t ) < 1.8 x 10~ 4 . We 
calculate the product of the B meson decay constant and \V u b\ to be • | | = (7.0^3g(stat.)^g 4 (syst.)) x 
10" 4 GeV. 
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1 INTRODUCTION 



In the Standard Model (SM), the purely leptonic decay B + —> t + v t 5 proceeds via quark annihilation into 
a W + boson (Fig. 1). Its amplitude is thus proportional to the product of the S-decay constant fs and the 
quark-mixing-matrix element V u b- The branching fraction is given by: 

2 

W1IK&I 2 , (i) 

where we have set % = c = 1, Gf is the Fermi constant, V u b is a quark mixing matrix element [1, 2], 
fs is the B + meson decay constant which describes the overlap of the quark wave-functions inside the 
meson, t b + is the B + lifetime, and m B and m T are the B + meson and r masses. This expression is 
entirely analogous to that for pion decay. Physics beyond the SM, such as a two-Higgs doublet models, 
could enhance or suppress the B{B + — > t + v t ) through the introduction of a charged Higgs boson [3]. 

Current theoretical values for (obtained from lattice QCD calculations) [4] have large uncertainties, 
and purely leptonic decays of the B + meson may be the only clean experimental method of measuring fs 
precisely. Given measurements of \V u b\ from semileptonic B — > ulv decays, fs could be extracted from 
the measurement of the B + — ► t + v t branching fraction. In addition, by combining the branching fraction 
measurement with results from B mixing, the ratio |V^&|/|Vfci| can be extracted from B{B + — > t + v t )/ Am, 
where Am is the mass difference between the heavy and light neutral B meson states. 



B(B^ 



G F m B m 2 
8^ 



mz 



m. 




Figure 1: The purely leptonic B decay B + — > t + v t proceeding via quark annihilation into a W + boson. 

The decay amplitude is proportional to the lepton mass and as such decay to the lighter leptons is 
suppressed. This mode is the most promising for discovery of leptonic B decays. However, experimental 
challenges such as the large missing momentum from several neutrinos make the signature for B + — ► t + u t 
less distinctive than for other leptonic modes. 

The SM estimate of this branching fraction is (1.59±0.40) x 1CT 4 , using \V ub \ = (4.39±0.33) x 1(T 3 [5] 
and f B = 0.216 ± 0.022 GeV [4] in Eq. 1. 

In a previously published analysis of a smaller sample of 223 x 10 6 T(4S) decays the BABAR collabo- 
ration set an upper limit of: 

B(B + -► t + u t ) < 2.6 x 10~ 4 at the 90% CL. [6] (2) 

5 Charge-conjugate modes are implied throughout this paper. The signal B will always be denoted as a B + decay while the 
semi-leptonic B will be denoted as a B~ to avoid confusion. 
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The Belle collaboration reported evidence of a signal in this channel recently [7]; a central value of 
(1.06lQ2g(stat.)^Q'}g(syst.)) x 10 -4 was extracted. The analysis presented herein is a detailed update of 
the previous BABAR search. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. The sample 
corresponds to an integrated luminosity of 288 fb -1 at the T(AS) resonance (on-resonance) and 27.5 fb _1 
taken 40 MeV below BB threshold (off -resonance). The on-resonance sample consists of about 320 x 10 6 
T(AS) decays (BB pairs). The collider is operated with asymmetric beam energies, producing a boost of 
/?7 0.56 of the T(45) along the collision axis. 

The BABAR detector is optimized for asymmetric energy collisions at a center-of-mass (CM) energy 
corresponding to the T(AS) resonance. The detector is described in detail in Ref. [8]. The components 
used in this analysis are the tracking system composed of a five-layer silicon vertex detector and a 40-layer 
drift chamber (DCH), the Cherenkov detector (DIRC) for charged ir-K discrimination, the Csl calorimeter 
(EMC) for photon and electron identification, and the 18-layer flux return (IFR) located outside of the 1.5T 
solenoidal coil and instrumented with resistive plate chambers for muon and neutral hadron identification. 
For the most recent 51 fb _1 of data, a portion of the muon system has been upgraded to limited streamer 
tubes (LST) [9]. We separate the treatment of the data to account for varying accelerator and detector 
conditions. "Runs 1-3" corresponds to the first 111.9 fb _1 , "Run 4" the following 99.7ftr 1 and "Run 5" 
the subsequent 76.8 fb _1 . 

A GEANT4-based [10] Monte Carlo (MC) simulation is used to model the signal efficiency and the 
physics backgrounds. Simulation samples equivalent to approximately three times the accumulated data 
were used to model BB events, and samples equivalent to approximately 1.5 times the accumulated data 
were used to model continuum events where e + e~ — > uu, dd, ss, cc and t + t~. A large sample of signal 
events is simulated, where a B + meson decays to t + v t and a B meson decays to an acceptable B mode. 
Beam related background and detector noise from data are overlayed on the simulated events. 

3 ANALYSIS METHOD 

Due to the presence of multiple neutrinos, the B + — > t + u t decay mode lacks the kinematic constraints 
which are usually exploited in B decay searches in order to reject both continuum and BB backgrounds. 
The strategy adopted for this analysis is to reconstruct exclusively the decay of one of the B mesons in the 
event, referred to as "tag" B. The remaining particle(s) in the event, referred to as the "signal side", are 
then compared with the signature expected for B + — > t + u t . In order to avoid experimenter bias, the signal 
region in data is not examined ("blinded") until the final yield extraction is performed. 

The tag B is reconstructed in the set of semileptonic B decay modes B~ — > D°£~i>£X, where £ is 
e or n and X can be either nothing or a transition particle from a higher mass charm state decay which 
we do not attempt to reconstruct (although those tags consistent with neutral B decays are vetoed). The 
D° is reconstructed in four decay modes: K~ir + , K~tt + it~it + , K~Tr + ir°, and i\~S?7r + 7r - . The K® is 
reconstructed only in the mode K® — > it + it~. These cases where the low momentum transition daughter 
of D*° decays need not be reconstructed and the final state B — ► D°£uX as observed provides a higher 
efficiency but somewhat lower purity than the exclusive reconstruction method of B~ — ► D £~vg. The 
choice of reconstructing the tag B as B~ — > D°£~i>iX was optimized by maximizing s/y/s + b where 
s = signal and b = background where a branching fraction for B + — > t + u t of 1 x 10~ 4 is assumed. 



9 



The B + — ► r + i/ r signal is searched for in both leptonic and hadronic r decay modes: r + — > e + f e i/ T , 
t + — ► ii + u^u T , t + — > 7r + i/ r and r + — ► tt + tt°u t . The branching fractions of the above r decay modes are 
listed in Table 1. 

Table 1: Branching fractions for the r decay modes used in the B + — > r + ^ T search [11]. 



Decay Mode 


Branching Fraction (%) 


r + — > e + v e v T 


17.84 ± 0.06 


r+ -»■ (l + V,j,U T 


17.36 ± 0.06 


T + — ► 7T + P r 


11.06 ±0.11 


T + — ► 7T + 7T l/ T 


25.42 ±0.14 



3.1 Tag B Reconstruction 

The tag 5 reconstruction proceeds as follows. First we reconstruct the D° candidates in the aforementioned 
four decay modes using reconstructed tracks and photons where a tt° is included. The tracks are required 
to meet particle identification criteria consistent with the particle hypothesis, and are required to converge 
at a common vertex. The 7r° candidate is required to have invariant mass between 0.115-0.150 GeV/c 2 and 
its daughter photon candidates must have a minimum energy of 30 MeV. The mass of the reconstructed 
D° candidates in K~ir + , K~ir + ir~ir + , and K^^tt^ modes are required to be within 20 MeV/c 2 of the 
nominal mass [11]. In the K~ir + ir° decay mode the mass is required to be within 35 MeV/c 2 of the nominal 
mass [11]. 

Finally D°£ candidates are reconstructed by combining the D° with an identified electron or muon with 
momentum above 0.8 GeV/c in the CM frame. The D° and £ candidates are required to meet at a common 
vertex. An additional kinematic constraint is imposed on the reconstructed D°£ candidates: assuming that 
the massless neutrino is the only missing particle, we calculate the cosine of the angle between the D°£ 
candidate and the B meson, 

2E B E D o e - m\ - m 2 D0e 

cos9 B _ DH = — ; (3) 

2\PB\\PD°e\ 

Here (E D t>£, Pc-op) an d (Eb, Pb) are the four-momenta in the CM frame, and m D 0£ and tub are the masses 
of the D°£ candidate and B meson, respectively. Eb and the magnitude of pb are calculated from the 



beam energy: Eb = -Ecm/2 and \ps\ = y E B — m? B , where Eb is the B meson energy in the CM frame. 
Correctly reconstructed candidates populate the range [—1,1], whereas combinatorial backgrounds can take 
unphysical values outside this range. We retain events in the interval —2.0 < cosO B _ D o£ < 1.1, where the 
upper bound takes into account the detector resolution and the loosened lower bound accepts those events 
where a soft transition particle from a higher mass charm state is missing. 

If more than one suitable D°£ candidate is reconstructed in an event, the best candidate is taken to be the 
one with the largest vertex probability. The sum of the charges of all the particles in the event (net charge) 
must be equal to zero. 

At this stage of the selection, the observed yield in data and the predicted yield in the MC simulation 
agree to within approximately 3%. This discrepancy is corrected by scaling the yield and efficiency obtained 
from MC simulation. By multiplying the relevant branching fractions and reconstruction efficiencies, from 
signal MC simulation, B tagging efficiencies are extracted. Scale factors of 1.05, 1.00 and 0.97 are used 
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to correct these efficiencies for Runs 1-3, Run 4 and Run 5 respectively. The systematic error associated 
with this correction is described in Sec. 5. The corrected tag reconstruction efficiency in the signal MC 
simulation is (7.61 ± 0.05) xl(T 3 for Runs 1-3, (6.31 ± 0.05) xl0~ 3 for Run 4 and (5.87 ± 0.06) xl0~ 3 
for Run 5 where the errors are statistical only. 

3.2 Selection of B + — > r + u T signal candidates 

After the tag B reconstruction, in the signal side the r from the B + — > t + v t decay is identified in one of 
the following modes: r + — ► e + u e u T , t + — > fi + u ll u T , t + — > 7r + i> r or r + — > 7T + 7r°i/ T . We select events with 
one signal-side track which must satisfy the following selection criteria: it must have at least 12 DCH hits, 
its momentum transverse to the beam axis, pt, is greater than 0.1 GeV/c, and its point of closest approach 
to the interaction point is less than 5.0 cm along the beam axis and less than 1.5 cm transverse to the beam 
axis. The invariant mass of a signal-side 7r° candidate must be between 0.115-0.150 GeV/c 2 , the shower 
shape of the daughter photon candidates must be consistent with an electromagnetic shower shape and the 
photons must have a minimum energy of 50 MeV in the CM frame. 

The different signal tau decay modes are distinguished by their selection criteria. The r + — > e + v e v T , 
r + — ► h+v^Vt, t + — > tt + u t and r + — ► 7r + 7r°i/ r signal modes, all of which contain one charged track, are 
separated by particle identification. Both the r + — > ir + i> T and the r + — > tt + tt°v t modes contain a pion 
signal track and are characterized by the number of signal-side ir° mesons. 

• Particle identification: 

- For the r + — > e + v e u T selection the track must be identified as an electron and not identified as 
a muon. 

- For the r + — > ii + u^i> T selection the track must be identified as a muon and not identified as an 
electron. 

- For the r + — ► it + v T selection we require that the track is not identified as an electron or a muon. 

- For the r + — > 7r + 7r°P T selection we require that the track is not identified as an electron or a 
muon or a kaon. 

• Signal-side 7r° multiplicity: 

- For the r + — > tt + u t selection we require the event to contain no signal-side tt°. 

- For the r + — > ir + ir v T selection we require that the event contains at least one signal-side ir°. 

Background consists primarily of B + B~ events in which the tag B meson has been correctly reconstructed 
and the recoil side contains one signal candidate track and additional particles which are not reconstructed 
by the tracking detectors or calorimeters. Typically these events contain K\ candidates and/or neutrinos, 
and frequently also additional charged or neutral particles which pass outside of the tracking and calorime- 
ter acceptance. Background events also contain B°B° events. The continuum background contributes to 
hadronic r decay modes. In addition some excess events in data, most likely from two-photon and QED 
processes which are not modeled in the MC simulation, are also seen. These backgrounds have a distinc- 
tive event shape and are suppressed by the following constraints on the kinematics of the B + — > t + u t 
candidates. 
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Missing mass: The missing mass is calculated as follows. 



M miss = yj (E T{4S) - E vis ) 2 - (p r( 4 5) - p vis ) 2 . (4) 

Here (E r ^ AS y Pr(4S)) i s tne four-momentum of the T(4S), known from the beam energies. The 
quantities E V - 1S and p v - ls are the total visible energy and momentum of the event which are calculated 
by adding the energy and momenta, respectively, of all the reconstructed charged tracks and photons 
in the event. 

> e + u e u T selection events with missing mass between 4.6 and 6.7 GeV/c 2 are 

► // + i/ jU i/ T selection events with missing mass between 3.2 and 6.1 GeV/c 2 are 

ir + & T selection the missing mass is required to be greater than 1.6 GeV/c 2 . 
7r + 7r°i/ T selection the missing mass is required to be less than 4.6 GeV/c 2 . 

Maximum CM momentum of the r daughter: 

The following maximum CM momentum requirements are applied to the r daughter particles. 

- The electron candidate from the r + — ► e + v e v T decay must have a CM momentum of less than 
1.5 GeV/c. The CM momentum requirement is not applied to the r + — ► n + v^v T selection 
because the momentum spectrum of the muon from r decays peaks below 1 GeV/c and the 
particle identification efficiency for low momentum muons is lower than that for low momentum 
electrons. Therefore, applying the maximum momentum cut reduces the selection efficiency of 
the r + — > /j, + i> l j i u T mode significantly. 

- For the two hadronic r decay modes, the CM momentum of the ir from r + — > tt + u t must be 
greater than 1.6 GeV/c. The tttt combination from r + — > 7r + 7r z/ T must have CM momentum 
greater than 1 .7 GeV/c. 

Continuum Rejection using the R TT variable: 

An effective way to remove e + e _ — > r + r~ background is to place a cut in a plane defined by two 
variables: the cosine of the angle between the signal candidate and the tag B's thrust vector (in the CM 
frame), and the minimum invariant mass constructable from any three tracks in an event (regardless 
of whether they are already used in a tag or signal candidates). For the background, the cosine of 
the thrust angle peaks at —1 and 1, while the minimum invariant mass peaks below 1.5 GeV/c 2 . We 
transformed this 2-D variable into a 1-D variable using the following empirically derived equation 

Rrr = y/(3.7 - \ co8(^ i8ignal ) | f + (Mf a - 0.75) 2 , (5) 

where M™" 1 is the minimum invariant mass of any three changed tracks and Of signal ^ s tne an §l e 
between the thrust axes of the reconstructed D£ and the signal candidates. Because other continuuum 
backgrounds also peak in the cosine of the thrust angle, this variable is good at rejecting other similar 
categories of non-66 background. The selection criteria imposed on this quantity are: 

- For r + -f e + v e v T : 2.78 < R TT < 4.0 



- For the r 
selected. 

- For the r 
selected. 

- For the r H 

- For the -H 
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- For t+ -»■ p + v p v T : R TT > 2.74 

- For r+ -»• vr+zV i? rr > 2.84 

- For r+ -► 7r + 7r°iv i? rr > 2.94 



The r + — > 7r + 7r°i/ r decay proceeds via an intermediate resonance. For this mode further background 
rejection can be achieved by applying the following requirements on the intermediate meson. 

- p + selection: 

The signal-side track is combined with a signal-side 7r° to form the p + candidate. In events with 
more than one signal-side 7r°, the candidate with invariant mass closest to the nominal tt° mass 
[11] is chosen. The invariant mass of the reconstructed p + is required to be within 0.64-0.86 
GeV/c 2 . A quantity similar to cos 9 B _ D og, which is defined in section 3.1, can be reconstructed 
for t — ► pv as follows: 

2E T E P -ml -ml 

cos 9 T . P = p - , (6) 

2\p 7 - 



where (E T , p T ) and (E p , p p ) are the four-momenta in the CM frame, m T and m p are the masses 
of the r and p candidate, respectively. The quantities \p T | and E T are calculated assuming the 
r is from the B + — > t + v t decay, and the B + is almost at rest in the CM frame. We accept 
candidates with cos6 T ^ p > 0.87. 

• -Eextra requirement: 

The most powerful variable for separating signal and background is the remaining energy (E extra ), 
calculated by adding the CM energy of the neutral clusters and charged tracks that are not associated 
with either the tag B or the signal. The photon candidates contributing to the E cxtra variable have 
minimum cluster energies of 20 MeV in the CM frame. For signal events the neutral clusters con- 
tributing to -Eextra arise predominantly from processes such as beam-background, hadronic split-offs 
and Bremsstrahlung. Signal events tend to peak at low E cxtrai values whereas background events, 
which contain additional sources of neutral clusters, are distributed towards higher -E ex tra values. The 
most signal sensitive region is optimized for each mode and is blinded in on-resonance data until the 
selection is finalized. The -E cx tra < 0.5 GeV region is defined as the nominal blinding region which 
is slightly larger than the signal region for each mode. 

For all the signal modes -E cx tra is optimized for the best signal significance (assuming the branching 
fraction is 1 x 10~ 4 ). The optimization yields to following requirements: 

- For r + — > e + v e v T : E cxtra <0.31 GeV 

- For r+ -»■ p + v p v T '- -Eextra <0.26 GeV 

- For r+ -► 7r+P r : E cxtra <0.48 GeV 

- For T+ -► TT + TT°U T : Eextra <0.25 GeV 

The signal selection criteria for all signal modes are summarized in Table 2. 
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Table 2: The selection criteria for different signal modes using a B — > D°£ tag are listed in this table. 











4.6 < M mlss < 6.7 


3.2 < M mlss < 6.1 


1.6 < M mlss 


M miss < 4.6 


Psiqnal ^ 1-5 




1-6 < Psiqnal 


^ Psiqnal 


No IFR 


2.78 < R TT < 4.0 


2.74 < R TT 


2.84 < i? TT 


2.94 < i? TT 


m ee > 0.1 GeV/ 








jV e if ra < 2 


N extra < 3 


N EMCK° < 2 










p* selection: 
0.64 < M p ± < 0.86 GeV 
0.87 < cos# T _p 


£ cxtra < 0.31 GeV 


£ cxtra < 0.26 GeV 


S cxtra < 0.48 GeV 


£ cxtra < 0.25 GeV 



Table 3: The signal efficiencies, mode-by-mode, relative to the number of tags. The branching fraction for 
the given r decay mode selected is included in the efficiency. 



Mode 


Efficiency (BF Included) 


t + — > e + v e 9 T 

T + — > 7T + P T 

r + — > 7r + 7r°P r 


0.0414 ± 0.0009 
0.0242 ± 0.0007 
0.0492 ± 0.0010 
0.0124 ± 0.0005 



3.2.1 Signal Efficiency 

The signal-side selection efficiencies for the r decay modes are determined from signal MC simulation and 
summarized in Table 3. The signal efficiencies correspond to the number of events selected in a specific 
signal decay mode, given that a tag B has been reconstructed. 

The selection efficiency for r + — > fi + v^i> T is low compared to that of the r + — > e + u e u T mode because 
the momentum spectrum of the signal muons peaks below 1 GeV/c, where the muon detection efficiency is 
low. Since no minimum momentum requirement and no tight pion identification criteria are applied to the 
r + — ► tt + u t signal selection, electron and muon signal tracks that fail particle identification requirement 
get selected in this mode. Any true r + — > ir + 7r u T signal events, with a missed 7r° also get included in 
t + — > 7t + pv selection mode. Therefore the r + — > 7r + z/ T selection mode has the highest signal efficiency. 

3.3 Validation of Background Estimation from -E ex tra Sidebands 

We further study the agreement between simulation and data by using the extra energy sideband region, and 
the ratio of the yields in this region to that in the signal region. This is used mainly to test the reliability of 
the background estimation in the low E cxtvll region by extrapolation from the higher E'extra region. 

The -Ecxtra > 0.5 GeV region is defined as the "sideband" (sb). The "signal region" is defined separately 
for each selection mode. For each control sample after applying appropriate selection cuts, the number 
of MC events in the signal region (iVMc.sig) and sideband (iVjvicsb) are counted and their ratio (R MC ) is 
obtained. 

-V M C,sb 
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Using the number of data events in the sideband (iVdaia,sb) and the ratio R MC , the number of expected 
background events in the signal region in data (-/Vgxp.sig) is estimated. 



^exp.Sig — -Wdata,sb " R (8) 

The number of expected data events (-/V exP; sig) in the signal region is compared with the observed number 
of data events (iVobs,sig) m the signal region. The agreement between the above two quantities provide 
validation of background estimation in the low E cxtJ -a region. 

Table 4 illustrates the level of agreement between the sideband projections in MC and data. In general, 
the agreement is at the la level between the direct count in the MC signal region and the projected data. 
The projections in data are used to predict background for the final extraction, hence we only rely on the 
data for this. 



Table 4: The sideband-to-signalbox projection computed using a sideband region where E cxtrau > 0.5 GeV. 
The second column corresponds to the ratio of yields in the signal region and sideband as measured in MC. 



Mode 


ratio (MC) 


upper sb (Data) 


signal region (Proj) 


signal region (MC) 


electron 


0.137 ± 0.015 


305.00 ±17.46 


41.91 ± 5.19 


39.72 ± 4.07 


muon 


0.037 ± 0.004 


965.00 ±31.06 


35.39 ±4.16 


36.13 ±4.02 


pion 


0.043 ± 0.004 


2288.00 ± 47.83 


99.09 ±9.10 


87.69 ± 7.72 


rho 


0.005 ± 0.001 


2805.00 ± 52.96 


15.30 ±3.48 


15.81 ± 3.58 



4 VALIDATION OF TAG B YIELD AND £ e xtra SIMULATION 

The tag B yield and E cxtTS , distribution in signal and background MC simulation are validated using various 
control samples. The level of agreement between the data and simulation distributions provides validation 
of the E'cxtra modeling in the simulation and corrects for differences in the yield of reconstructed tag B's. 

"Double-tagged" events, for which both of the B mesons are reconstructed in tagging modes, B~ — > 
D°l~viX vs. B + — ► D°l + vnX are used as the main control sample. Due to the large branching fraction and 
high tagging efficiency for these events, a sizable sample of such events is reconstructed in the on-resonance 
dataset. Due to all of the decay products of the T{AS) being correctly accounted for the double-tagged 
events reconstructed have a high purity. 

To select double-tag events we require that the two tag B candidates do not share any tracks or neutrals. 
If there are more than two such non-overlapping tag B candidates in the event then the best candidates 
are selected as those with the largest D°-£ vertex probability, as with the signal search. The number of 
double-tagged events (AT 2 ) is given by 

N 2 = e 2 N. (9) 

where N is the number of BB events in the sample and e is the tag efficiency that is compared between data 
and MC. Using the expression in equation 9 we calculate the efficiencies £data and £mc- The correction fac- 
tor, ratio of the efficiencies between data and simulation, from this method is given in equations 10, 1 1 and 12 
for Runs 1-3, Run 4 and Run 5 respectively. 

gRuns 1-3 = 1.05 ± 0.02 (10) 

£mc 
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=1.00 ±0.03 (11) 



=0.97 ±0.03 (12) 

It was directly verified that data taken during Runs 1-3 agreed in both shape and normalized yield whereas 
during Run 4 and Run 5 data were taken with the machine operating in a mode of continuous injection 
which may affect detector backgrounds differently. These runs are therefore considered separately. 

The E'oxtra for the double-tagged sample is calculated by summing the CM energy of the photons which 
are not associated with either of the tag B candidates. The sources of neutrals contributing to the -Eextra 
distribution in double-tagged events are similar to those contributing to the E extI3i distribution in the signal 
MC simulation. Therefore the agreement of the E extIa distribution between data and MC simulation for the 
double-tagged sample, in figure 2, is used as a validation of the £7 e xtra simulation in the signal MC. 
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Figure 2: The distribution of the remaining neutral energy (E extra ) for double-tagged events, plotted for 
generic MC and data: a) Runs 1-3, b) Run 4 and c) Run 5. No off-resonance data events are seen in the 
-Eextra region plotted here. In these events both of the D°£ candidates from double-tag are required to pass 
the selection described in section 3.1 and best candidate selection. The differences in these distributions are 
used for obtaining the systematic error for tagging efficiency correction. 

The simulation is further validated by comparing a sample of events where the signal candidate and tag 
B candidate are of the "wrong-sign" with non-zero net charge. The agreement between data and simulation 
for all signal modes for the background estimation in the E cx tra signal region provides a useful cross-check. 



5 STUDIES OF SYSTEMATICS 

The main sources of uncertainty in the determination of the B 
ing: 



t + v t branching fraction are the follow - 



• Uncertainty in tagging efficiency determination 

• Uncertainty in determination of the efficiency for each selection mode. 

• Uncertainty in the determination of the number of expected background events in the signal region 
for each selection mode. 



A small uncertainty of 1 . 1 % also enters the branching ratio limit calculation from the estimation of the 
number of B + B~ events present in the data sample [12]. The systematic uncertainties are summarized in 
table 5. 
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5.1 Tagging Efficiency Systematics 

The tagging efficiency and yield in signal simulation is corrected using the double-tagged events. The 
selection of double-tagged events is described in section 4. 

We take the 1.9%, 3.0% and 3.1% errors (from equations 10, 11 and 12) obtained from the double tag 
method as the systematic uncertainties associated with the tagging efficiency and yield correction in MC. 
The combined, luminosity weighted, tag B yield systematic uncertainty is 1.5%. The luminosity weighted 
tag B yield correction is 1.01. 

5.2 E extra Systematic Uncertainty 

The systematic uncertainty due to the mis-modeling of the E cx tra variable is extracted using the double- 
tagged events. The selection of double-tagged events is described in Section 4. A cut is imposed on 
the -Eextra distributions shown in Figures 2(a), 2(b) and 2(c) to extract the yield of candidates satisfying 
Eextra < 0.5GeV. This yield is then compared to the number of candidates in the full sample. Comparing 
the ratio extracted from MC to that extracted from data yields a correction factor, the error on which is 
taken as the systematic uncertainty for E ex tra- These values are broken up by run and we extract the follow- 
ing numbers: Runs 1-3 = 0.98±0.06, Run 4 = 0.99±0.06, Run 5 = 1.02±0.08 The combined, luminosity 
weighted systematic uncertainty for E ex tra is 3.8%. The luminosity weighted E cx tra correction is 0.99. 

5.3 Uncertainties in the signal selection efficiencies in each selection mode 

Besides the tagging efficiency uncertainty, the contribution to the systematic uncertainties in the determina- 
tion of the efficiencies comes from systematic uncertainty on the tracking efficiency, particle identification, 
and simulation of the neutral clusters in the calorimeter which contribute to the Eextra distribution, and 
K L identification. The different contributions to the systematic uncertainty on the selection efficiencies are 
listed in table 5. 

Table 5: Contribution to the systematic uncertainty on the signal selection efficiencies in different selec- 
tion modes. These uncertainties are added together in quadrature with the uncertainty on the tag B yield, 
extracted from the double-tagged control sample, of 1.5%. The uncertainty on MC statistic is added in 
quadrature to obtain the total systematic uncertainty. 



Selection 


tracking 


Particle 


K 


E e xtra 


7T U 


Total 


Correction 


modes 


(%) 


Identification 




modeling 


modeling 


Systematic 


Factor 






(%) 


(%) 


(%) 


(%) 


Error (%) 




e r v e y T 


0.3 


2.0 


3.6 


3.8 




5.8 


0.982 




0.3 


3.0 


3.6 


3.8 




6.2 


0.893 


TT + U T 


0.3 


1.0 


6.2 


3.8 




7.5 


0.966 




0.3 


1.0 


3.6 


3.8 


1.8 


5.8 


0.961 



5.4 Uncertainties on K° L modeling 

The systematic uncertainty on the modeling of K® candidates is extracted using the double-tagged events 
outlined in section 4. A comparison between data and simulation is used to extract both a correction and 
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a systematic uncertainty, similarly to the method used for -E^tra- We quantify this comparison by com- 
paring the yield with a cut demanding exactly zero reconstructed IFR measured K\ candidates remaining, 
with a sample where any number of K ° candidates remain and take the ratio of ratios from the MC and 
data. We extract the following values for corrections and systematic uncertainties: Runs 1-3 = 0.98±0.05, 
Run 4 = 1.00±0.06, Run 5 = 0.98±0.08, hence percentage uncertainties of 5.1%, 6.0% and 8.2%. The cor- 
rection factors are all close to unity as expected. The combined, luminosity weighted systematic uncertainty 
for IFR K® candidates is 3.6%. The luminosity weighted IFR K® correction is 0.99. 

The same exercise is performed for K® candidates reconstructed in the EMC. We extract the fol- 
lowing values for corrections and systematic uncertainties: Runs 1-3 = 0.88±0.05, Run 4 = 1.00±0.10, 
Run 5 = 1.08±0. 11. Percentage uncertainties are 5.7%, 10% and 10.2%. The combined, luminosity 
weighted systematic uncertainty for EMC K® candidates is 5.1%. The luminosity weighted EMC K® 
correction is 0.97. 



6 RESULTS 

After finalizing the signal selection criteria, the signal region in the on-resonance data is examined. Table 6 
lists the number of observed events in on-resonance data in the signal region, together with the expected 
number of background events in the signal region. Figures 3 and 4 show the E'extra distribution in data 
and simulation for each of the r decay modes considered. Data is overlayed on the summed MC contribu- 
tion, scaled to the dataset luminosity, and signal MC is plotted for comparison. Figure 5 shows the -E ex tra 
distribution for all modes combined. 
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Figure 3: Total extra energy is plotted after all cuts have been applied in the mode (a) r + — > e + u e u T and 
(b) r + — ► fi + u^i' T . Off -resonance data and MC have been normalized to the on-resonance luminosity. 
Simulated B + —* t + v t signal MC is plotted (lower) for comparison. 
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Figure 4: Total extra energy is plotted after all cuts have been applied in the mode (a) r + — > ir + v T and (b) 
r + — ► 7r + 7r°z/ T . Off-resonance data and MC have been normalized to the on-resonance luminosity. Simu- 
lated B + — ► t + u t signal MC is plotted (lower) for comparison. 



We determine the B + — » t + v t branching fraction from the number of signal candidates S{ in data for 
each r decay mode, according to S{ = N B -gB{B + — * r + i/ T )e tag ej. Here N B -g is the total number of BB 
pairs in data, e tag is the tag reconstruction efficiency in signal MC; e% is the signal-side selection efficiency in 
different r decay modes calculated with respect to the total number of reconstructed tag B mesons. Table 7 
shows the values of N B g, Stag and e- L after applying appropriate systematic corrections (see section 5). The 
results from each decay mode are combined using the ratio Q = C(s + b)/C(b), where C(s + b) and C(b) 
are the likelihood functions for signal plus background and background-only hypotheses, respectively [13]: 

C(s + b) s n m nn^- d3) 

1=1 % 1=1 1 

We include the statistical and systematic uncertainties on the expected background (bi) in the likelihood 
definition by convolving it with a Gaussian distribution {Q). The mean of Q is bi, and the standard deviation 
((Jbi) of Q is tne statistical and systematic errors on bi added in quadrature [14], 

£{si + bi) ^ £( Si + bi)®g(bi,a bi ) (14) 

(similarly for C(b{)). The results from this procedure are illustrated in Figure 6. 
We determine the following branching fraction 

B(B + -» t + v t ) = (0.88lg;g?(stat) ± O.ll(syst)) x 10~ 4 , (15) 

and also set an upper limit at the 90% confidence level of 

B(B + -> t + u t ) < 1.8 x 10" 4 . (16) 
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Figure 5: Total extra energy is plotted after all cuts have been applied with all modes combined. Off- 
resonance data and MC have been normalized to the on-resonance luminosity. Events in this distribution 
are required to pass all selection criteria. In addition the background MC have been scaled according to the 
ratio of predicted backgrounds from data and MC as presented in section 3.3. Simulated B + — > t + v t signal 
MC is plotted (lower) for comparison. 



Figure 6 shows the distributions of confidence level vs branching fraction and the negative log likelihood 
curve illustrating the extracted upper limit and central value respectively. 

Using the measured central value for B(B + — ► t + u t ) and taking the known values of Gf, tub, m T and tb 
from Ref. [11] we calculate, from equation 1, the product of the B meson decay constant and \ V u b\ to be 
Sb ■ \V ub \ = (7.0l!;jj(staL)18;f (syst)) x 10~ 4 GeV. 



7 SUMMARY 

We have performed a search for the decay process B + — > t + v t . To accomplish this a sample of semileptonic 
B decays (D°£~i>£X) has been used to reconstruct one of the B mesons and the remaining information in 
the event is searched for evidence of B + — > t + u t . A branching fraction of 

B{B + -► t + v t ) = (0.88lg;P(stat.) ± O.ll(syst)) x 10~ 4 , (17) 

is measured and we set an upper limit at the 90% confidence level of 

B(B + t + u t ) < 1.8 x 10~ 4 . (18) 

Using the measured central value for B(B + — > t + v t ) and taking the known values of Gf, tub, m T and tb 
from Ref. [11] we calculate, from equation 1, the product of the B meson decay constant and |V^{,| to be 
Sb ■ \V ub \ = (7.0l|| (staL)lg;f (syst)) x 10~ 4 GeV. 
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Table 6: The observed number of on-resonance data events in the signal region are shown, together with 
number of expected background events. The background estimations include systematic corrections referred 
to in section 3.3. 



Selection 


Expected 


Observed Events 




Background Events 


in On-resonance Data 


e + v e v T 


41.9 ± 5.2 


51 




35.4 ± 4.2 


36 


TT + U T 


99.1 ±9.1 


109 




15.3 ±3.5 


17 


All modes 


191.7 ± 11.8 


213 



Table 7: The corrected tag and signal efficiencies. Two errors are quoted: the first is the MC statistical 
uncertainty, and the second is the systematic error computed from the sources in section 5. 



Efficiency 


Corrected 


Relative Systematic Error (%) 


Tag 


(6.77 ± 0.05(stat.) ± O.lO(syst)) x 10" 3 


1.5 


e(r + — > e + v e v T ) 


(4.06 ± 0.09(stat.) ± 0.23(syst.)) x 10" 2 


5.6 


e(r+ -> ii+v^Vr) 


(2.16 ± 0.06(stat.) ± 0.13(syst.)) x 1Q-' 2 


6.0 


e{t + — > TT + V T ) 


(4.88 ± O.lO(stat) ± 0.35(syst.)) x 10" 2 


7.3 


E{T + — > 7T + 7T U i/ T ) 


(1.16 ± 0.05(stat.) ± 0.07(syst.)) x 10~ 2 


5.6 




Branching Fraction (10 4 ) Branching Fraction (10 -4 ) 



Figure 6: The confidence level vs branching fraction is shown (left) to illustrate the extracted upper limit. 
The negative log likelihood curve (right) illustrates the central value and it's corresponding uncertainty. 
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